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Abstract 

The accurate determination of the elemental composition of cosmic rays at high 
energies is expected to provide crucial clues on the origin of these particles. Previ- 
ous direct measurements of composition have been limited by experiment collecting 
power, resulting in marginal statistics above 10 14 eV, precisely the region where 
the "knee" of the cosmic-ray energy spectrum is starting to develop. In contrast, 
indirect measurements using extensive air showers can produce sufficient statistics 
in this region but generate elemental measurements which have relatively large un- 
certainties. Here we discuss a technique which has become possible through the use 
of modern ground-based Cerenkov imaging detectors. We combine a measurement 
of the Cerenkov light produced by the incoming cosmic-ray nucleus in the upper 
atmosphere with an estimate of the total nucleus energy produced by the exten- 
sive air shower initiated when the particle interacts deeper in the atmosphere. The 
emission regions prior to and after the first hadronic interaction can be separated 
by an imaging Cerenkov system with sufficient angular and temporal resolution. 
Monte Carlo simulations indicate an expected charge resolution of /S.Z/Z < 5% for 
incident iron nuclei in the region of the "knee" of the cosmic-ray energy spectrum. 
This technique also has the intriguing possibility to unambiguously discover nuclei 
heavier than iron at energies above 10 14 eV. The identification and rejection of 
background produced by charged particles in ground based gamma-ray telescopes 
is also discussed. 
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1 Introduction 



The origin of cosmic rays remains a central unresolved question in astrophysics. 
Nearly 90 years after the first observations, this population of charged particles 
remains an enigma; the heart of which is the huge dynamic range of fluxes 
and energies over which they have been observed. In the theoretical arena, 
a relatively recent paradigm has emerged, which involves the acceleration of 
cosmic rays via diffusive shock processes in supernovae remnants (SNR). This 
idea has been fueled both by a viable physical model of the acceleration process 
[1-4] and a simple energetics argument in which galactic supernovae are the 
only galactic candidate with sufficient energy output to supply the staggering 
amount of power needed to sustain the cosmic-ray population (see e.g., [5]). 

A key issue with the SNR idea is that supernova diffusive shock acceleration 
can only produce particles up to some maximum energy, limited either by the 
lifetime of the strong shock or by the particles becoming so energetic they can 
no longer be confined in the acceleration region [6]. Estimates of this upper 
energy limit vary, with typical values in the region of 10 14 eV. However, the 
observed flux of cosmic rays extends more or less continuously for another five 
orders of magnitude beyond this; no plausible extrapolation of the standard 
parameters of SNR can generate cosmic rays of such energies. 

Another aspect of the cosmic ray riddle is the existence of an observational 
"knee" at ~ 10 15 eV in the cosmic ray energy spectrum. The coincidence of 
this "knee" with the theoretical energy limit of SNR diffusive shock acceler- 
ation is intriguing and has actually served as evidence to support the theory. 
However the "knee" represents only a small change in the spectral slope of 
the overall flux, with the energy dependence changing from E~ 2 - 75 below the 
"knee" to E~~ 3 - above it. We are faced with a simple observational fact that 
the cosmic rays have an essentially continuous spectral slope for nearly 11 
orders of magnitude. 

In principle, additional mechanisms could provide the flux at high energies. 
The power budget for "post-knee" cosmic rays is only a fraction that of the 
total budget, so there is some freedom in selecting models. However, to result 
in an energy spectrum as smooth as is observed, these mechanisms would 
have to generate fluxes which are remarkably, perhaps implausibly, close to 
that of the SNR mechanism. Another problem exists at the highest energies, 
around ~ 10 19 eV, where an additional spectral break occurs (the "ankle"). 
It has been argued that particles above the "ankle" could be extragalactic. If 
so, they can be expected to exhibit a characteristic cutoff due to photo-pion 
production with the cosmic microwave background. This cutoff has yet to be 
observed [7]. 
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An overall scheme which credibly addresses and unifies these issues remains 
elusive. More reliable and accurate measurements in the "knee" region, of 
abundant cosmic ray nuclei (1 < Z < 26) could drastically revise our current 
ideas and provide a basis for such a unification. The accurate determination of 
the composition of cosmic rays has provided some of the key advances in this 
field at lower energies, where direct measurements are possible with detectors 
above the atmosphere. For example, the realization that the observed spectral 
slope of cosmic rays is significantly steeper than that produced in the cosmic 
ray sources themselves resulted from measurements with sufficient elemental 
resolution to separate primary source cosmic ray elements from those produced 
in the interstellar medium at 10 11 eV [8]. 

In this paper we propose an idea which has become possible through advances 
in the imaging atmospheric Cerenkov technique. It is based on the concept that 
a detector of fine enough pixelation will be capable of observing Cerenkov 
light emitted directly from cosmic ray nuclei prior to their first interaction 
in the atmosphere. In general, this light is overwhelmed by the Cerenkov 
emission from the subsequent extensive air shower (EAS). However, with an 
appropriate detector, and within certain geometric constraints, this "direct 
Cerenkov" , or DC, light can be sufficiently well separated from the background 
of EAS Cerenkov light to make relatively high-precision measurements of the 
traditional cosmic ray composition (1 < Z < 26). This technique may also 
provide new opportunities for the first measurement of higher charge nuclei (Z 
S> 26) as well as possible improvements in ground based gamma-ray telescope 
sensitivity. 



2 Method 

2. 1 History 

Cerenkov light from extensive air showers was first predicted by Blackett in 
1948 [9], and later observed by Galbraith and Jelley in 1952 [10]. In 1989, the 
Whipple telescope used the technique of Cerenkov light imaging to provide the 
first highly significant detection of the Crab Nebula in high energy gamma rays 
[11]. The imaging technique, which was proposed in 1977 [12], involves the use 
of an array of photomultiplier tubes at the focal plane of the telescope to reject 
the very large background of cosmic ray-induced air showers by the shape of 
their images in the field of view. A camera with fine enough pixelation has 
the capacity to not only discriminate between hadronic and electromagnetic 
showers, but to map the development of the showers as they penetrate into 
the atmosphere. 
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In the present work we propose a means to measure the direct Cerenkov (DC) 
light produced by the incoming nucleus prior to its first interaction by using the 
imaging atmospheric Cerenkov techniques of VHE gamma ray astronomy. The 
targeting of DC light is not a new idea. In 1965, Sitte [13] had proposed that 
this radiation might be observed in high-altitude balloon-borne instruments. 
His idea, which was revisited by Gough in 1976 [14], was to place detectors 
at a height in the atmosphere above the mean interaction point of heavy 
primary cosmic rays and to look for the direct production of Cerenkov light. 
In the absence of the large light backgrounds due to emission from EAS, it was 
expected that accurate composition measurements could be made by analysis 
of the direct Cerenkov light yields. After a pioneering flight by Sood in 1981 
[15], the concept was unexploited until recent efforts by Seckel et al. in 1998 
[16]. The fundamental challenge in the new approach discussed here is the 
identification of the DC light against the much larger background of Cerenkov 
light produced in the associated EAS. 

2.2 Cerenkov Radiation 

A charged particle traveling in the atmosphere will produce Cerenkov radiation 
if it has a velocity greater than the local velocity of light. The threshold Lorentz 
factor, 7q, at which the radiation starts to be emitted is approximately: 

where 5 — n — 1, and n is the local index of refraction in the atmosphere. At 
sea level, 70 ~ 42, while at an altitude of some 50 km, 70 ~ 680. For an Iron 
nucleus, this is equivalent to energies of 2 TeV and 36 TeV, respectively. 

The rate of emission of Cerenkov light, Ng, increases rapidly with 7 above 
this threshold as: 

<* Z 2 <4 - ( 2 ) 

where Z is the particle charge. For 7 ^> 70 the amount of emitted light ap- 
proaches a saturation level which scales exclusively with the square of the 
particle charge and the local atmospheric density. The Lorentz threshold, 70, 
is a function of the local density, which, in turn, is a function of the altitude. 
Therefore, because the atmospheric density scales with height, once in satura- 
tion, the amount of Cerenkov emission from a particle is determined entirely 
by the charge of the particle, and the altitude of emission. The Cerenkov an- 
gle too, is determined by the local density (and thus, on the emission height), 
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Fig. 1. Schematic Representation of the Cerenkov emission from a cosmic ray pri- 
mary. 

and so the light pool from an emitting particle has a well-defined geometry 
determined entirely by the atmospheric density profile. 

In contrast, the Cerenkov emission from an extensive air shower has a yield 
almost linearly proportional to the energy of the primary particle, with a weak 
dependence on the primary charge, and a geometry which varies primarily with 
the details of the shower development. Figure 1 shows a schematic depiction 
of the Cerenkov radiation from a primary cosmic ray and its subsequent EAS. 
In the interior regions of the Cerenkov light pool, the typical photon densities 
due to the EAS emission are many times greater than those associated with 
the direct emission. Thus, to discern the DC photons against this background, 
some distinguishing characteristics must be identified. 

2.2.1 Angular Characteristics 

One of the fundamental characteristics of atmospheric DC light is a simple 
correspondence between emission height and emission angle. As shown above, 
this stems from the density dependence of the Cerenkov emission angle and 
the height dependence of the atmospheric density. In practice, this means that 
for fixed geometries, single emitting particles can be matched by their emission 
angle to unique emission heights. 
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Figure 2 shows the characteristic features of DC light from a relativistic 
(7 ^> 70) charged particle vertically incident on the top of the atmosphere. To 
illuminate the features of the DC emission, the particle is constrained to not 
interact in the atmosphere; the intersection point of the particle trajectory 
with the observing plane is the point where the impact parameter is zero. The 
default CORSIKA Monte Carlo simulation atmospheric profile [17] is used. 

The upper panel shows the impact parameter of photons at sea level pro- 
duced as a function of the particle altitude in the atmosphere. The physics of 
atmospheric Cerenkov emission produces a narrow emission cone angle at high 
altitude which expands as the particle penetrates deeper into the atmosphere. 
At altitudes around 10 km, the interplay between increasing emission angle 
and decreasing emission height produce a pileup effect near the maximum 
geometrically-allowed radius, near ~ 145 m. This pileup is responsible for the 
so-called "Cerenkov ring" . 

Since the rate of emission of Cerenkov light increases with atmospheric depth 
(because 70 becomes smaller), the optimal impact parameter for viewing DC, 
in the absence of an EAS, would be at this Cerenkov ring. However, since 
most nuclei, in reality, interact by altitudes of 25 km or so, the best impact 
parameter must be chosen as a compromise between high DC density (low 
altitude), and low nuclei interaction probability (high altitude). A reasonable 
radius is ~ 80 m. At this radius, a 10 m diameter detector will view DC light 
coming from an angular range of ~ 0.15 degrees with respect to the incoming 
trajectory, which corresponds to heights between roughly 29 and 34 km in the 
atmosphere. On average, this is about as deep as the cosmic ray primary can 
be expected to penetrate without interaction. 

As Figure 2 shows, an additional component of secondary Cerenkov light which 
originates from emission at an altitude ~ 4 km could also be observed at 
80 m. However, this light is emitted at an angle of nearly 1 degree, which is 
distinguishable from the light produced at higher altitudes by using an imaging 
detector with fine enough pixels. In a telescope with 0.05 degree pixels, the 
DC light from high altitude would be confined to a specific range of pixels. 
In the absence of EAS particles, measurements of photon densities in these 
angular bins would provide an unambiguous determination of the DC yield 
from the primary particle, and therefore its charge. 

However, the presence of the EAS acts to obscure the direct Cerenkov com- 
ponent, as electrons produced at lower altitudes get scattered, producing 
Cerenkov light which enters into the angular bin of the direct radiation. In 
lower energy showers, the effect is small compared to the densities of the DC 
light. However, as the size of the EAS increases, the background grows to the 
point where the density of the scattered electron emission exceeds that of the 
DC emission. Because the amplitude of the EAS light is set by the energy of 
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Fig. 2. The direct Cerenkov emission characteristics of a single non-interacting par- 
ticle vertically incident on the atmosphere observed at sea level. Upper Panel: The 
Cerenkov photon emission height versus photon impact parameter relative to the 
original particle trajectory. Lower Panel: The photon time delay measured at the 
observation depth relative to the particle traveling at speed c versus photon impact 
parameter relative to the original particle trajectory 

the incoming nucleus, this background effectively provides an upper limit to 
the useful energy range for the DC technique. This energy limit essentially 
scales with the square of the primary nuclear charge. 



2.2.2 Temporal Characteristics 

The lower panel of Figure 2 shows the time delay of Cerenkov photons emitted 
by the same particle as seen in the first panel. The delays are shown as a func- 
tion of the photons' final impact parameters at sea level. These are measured 
with respect to the time it takes the particle (assumed to be traveling at c) to 
reach sea level. 
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The DC light has a very fast time structure at high altitude. All of the pho- 
tons emitted from altitudes higher than 30 km arrive within roughly 300 ps of 
each other. Furthermore, this light is delayed with respect to the light emitted 
lower in the atmosphere (e.g., below 5 km) by 3 ns or more, depending on the 
observation radius. The effect of optical dispersion over a typical wavelength 
range for atmospheric Cerenkov measurements (300 — 600 nm) has a negligible 
effect on these timing widths. This leads to a additional method for differen- 
tiating between the DC light and that produced in the EAS. The DC should 
be delayed from the main Cerenkov pulse, and should have a characteristic 
width over an order of magnitude shorter than the EAS Cerenkov pulse. The 
width of an EAS Cerenkov pulse depends on the width of the shower front. 
These fronts are typically rs 2 m thick, leading to EAS Cerenkov pulses on 
order 6 ns wide. 

The largest time separation between DC and EAS Cerenkov light occurs clos- 
est to the core of the shower. However the total amount of DC radiation 
collected is smaller closer to the core as discussed above. Again, an observa- 
tion radius of ~ 80 m is a good compromise between providing an adequate 
DC signal before the particle is likely to interact and maintaining some time 
separation between the DC and EAS signal as shown in the lower panel of 
Figure 2. 



3 Simulation 

Since the dominant background for the DC light is Cerenkov light from elec- 
trons scattered in the EAS development, we have used numerical simulations 
to study the levels and fluctuations in this background. The characteristics of 
DC light in EAS have been modeled using a modified MOCCA Monte Carlo 
Simulation [18-20] and also the CORSIKA (Version 5.945, QGSJet98) sim- 
ulation package [17]. Both of these codes demonstrate similar characteristics 
for the DC light, and predict that DC light should be observable against the 
background EAS Cerenkov light over an energy window which depends on the 
charge of the primary particle. 

3.1 Average Behavior 

First we consider the average behavior of the DC light signal from various 
incoming particles observed at sea level. The left hand panel of Figure 3 illus- 
trates the angular and time characteristics of DC light emitted by a single 10 
TeV vertical cosmic ray with Z — 12 (Mg) which interacts and produces an 
EAS deeper in the atmosphere. The time axis corresponds to the time delay 
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Fig. 3. Simulated DC and EAS Cerenkov light emitted from a single vertically 
incident particle. Left-hand panel shows a 10 TeV Z = 12 nucleus and the right 
hand panel shows a 5 TeV Z = 4 nucleus. The Cerenkov light is observed at a 
radius between 67 — 94 m (mean radius 80 m) from the shower axis. The vertical 
axis is the time delay of the arriving photons as discussed in the text. The horizontal 
axis is the arrival angle of the photons with respect to the vertical at the observing 
site. The scales below each panel give the photon intensities. 

compared to a reference arrival time at sea level of a particle traveling at a 
speed c along the incoming nucleus path. The angular axis corresponds to 
the angle of entry of the photons into the detector compared to the incoming 
particle trajectory. The photon density of the light is averaged over an annu- 
lus extending from 67 m to 94 m from the shower axis, giving an 80m mean 
radius of the annulus, and the intensity is integrated over a wavelength band 
of 300 — 600 nm. The DC light emission is clearly seen as an arc on the left of 
the figure, separated from the Cerenkov emission produced by the EAS. The 
right hand panel of Figure 3 shows the signal from a single 5 TeV Z = 4 (Be) 
nucleus. This shower also has an obvious DC signal, but the intensity of this 
feature is clearly less well separated from the background. 

Above a certain energy, the light produced by the EAS development will pro- 
vide a strong enough background to completely overwhelm the DC light in- 
tensity, thereby making reliable measurements impossible. The DC light can 
therefore only be determined over a limited energy window. The lower energy 
threshold is denned by the threshold for Cerenkov emission, set by 70 (Fig- 
ure 4), while the upper energy threshold is set as the energy at which the 
secondary EAS Cerenkov light overwhelms the DC light. 

For heavy nuclei the lower total particle energy threshold for the observation 
window, Ei, can be described by E\ = k±Z where k\ is a constant. The up- 
per energy limit for DC observation, E^, occurs when the DC light level is 
essentially equal to the background light level. Since DC light emission N D q 
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Fig. 4. Simulated Average DC Emission rate for different charge primary nuclei. All 
shower trajectories are vertical and measurements are made at sea level. Vertical 
Axis: Average Emission rate (Photons/m 2 at 80 m mean radius from shower core). 
Horizontal Axis: Primary cosmic ray Lorentz factor 7. 

is proportional to Z 2 , we can express this as N D( j = k 2 Z 2 . The EAS Cerenkov 
light Ns is proportional to E, A^ s = k^E. Combining these last two equations 
gives an expression for the upper energy observation limit: 



E h = (k 2 /k 3 )Z 2 

Hence the relative width of the observation window AE/ E\ = (Eh — E{)/ ' E\ is 



AE/Et = (h/hk^Z - 1 

This window expands like ~ Z for heavy nuclei. This arises because the DC 
light scales like Z 2 , whereas the background level scales like Z. The results 
of numerical simulations give similar results to this simple estimate. Figure 
5 shows the upper and lower energy limits for DC observation as a function 
of primary particle charge Z. A DC measurement of Iron nuclei (Z = 26) is 
possible into the cosmic-ray "knee" region, around 10 3 TeV. 

The separation between the DC light and the secondary Cerenkov light per- 
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Fig. 5. Upper and lower threshold for detection of DC light in cosmic ray air showers. 
The lower threshold is due to the Cerenkov photon emission threshold. The upper 
threshold is where secondary light density from the EAS has a strength equal to the 
DC light density. Vertical Axis: Primary Energy (TeV). Horizontal Axis: Primary 
particle charge Z. 

sists at relatively large cosmic ray zenith angles. The lower right panel of 
Figure 6 shows the image from a theoretical Z = 50 nucleus at 45° to the 
zenith. Since both the DC light radius and the size of the EAS Cerenkov disk 
grow proportionally with the distance to the observation level, the separation 
between these two light emission regions persists. An advantage in large angle 
observation is that the light pool size grows geometrically with zenith angle, 
thereby increasing the effective detection area for cosmic ray observation. This 
may prove important for increasing the detector collection area at the highest 
energies (large Z) where cosmic ray particle fluxes are expected to be low. 



4 Resolution Considerations 

Above threshold, the DC light emission is independent of primary energy, 
making it an ideal measure of the primary charge. However, due to various 
background, detector resolution, and fluctuation processes, what one actually 
measures is not only the DC light emission, but rather the integrated DC light 
yield plus all background across the length of the primary cosmic track. We 
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Fig. 6. Simulated DC and EAS Cerenkov light emitted from a single particle. The 
upper left hand panel shows a 10 TeV 7 at vertical incidence. The upper right panel 
is 100 GeV electron at vertical incidence. The lower left panel shows a 250 GeV 
proton at vertical incidence, and the lower right panel shows a 200 TeV Z = 50 
nucleus with a trajectory at 45° to the zenith. The Cerenkov light is observed at 
a radius between 67 — 94 m (mean radius 80 m) from the shower axis. The scales 
below each panel give the photon intensities. The vertical axis is the time delay of 
the arriving photons as discussed in the text. The horizontal axis is the arrival angle 
of the photons with respect to the vertical at the observing site. The scales below 
each panel give the photon intensities. 

describe each of these effects and their Z dependence, and present a composite 
plot of the predicted charge resolution effects as a function of the primary 
charge Z. 



o.o 



o 

w 
c 

10.0 




Q 

E 



20.0 - 




12 



3000 



2500 



2000 



1500 



1000 



500 



Z=26 E=100 TeV 



57 m 



80 m 



10 



20 



30 



40 



50 



60 



70 



80 



90 



100 



Photons/m 



Fig. 7. Histogram of DC photon yield observed at mean radii of 57 m and 80 m 
from shower core for a large number of 100 TeV iron nuclei (Z=26). Horizontal Axis: 
Photon density (photons/m 2 ). Vertical axis: Number of Events 



4-1 Background Light Levels 



Background light sources which may degrade the charge resolution of this 
technique include fluctuations in the night sky background level, fluctuations 
in the secondary Cerenkov light emitted by the extensive air shower, and 
secondary light which is scattered by the atmosphere to the same direction 
and arrival time as the DC light signal. 



4- 1.1 Night Sky Background 

The overall background expected by fluctuations in the night-sky light level 
can be determined from previous measurements at dark sites. The standard 
background value quoted is 2 x 10 12 photons/(m 2 sec sr) over the range 300 — 
600 nm. This translates to an overall rms intensity scale near 0.01 photons/m 2 
for the angle/time bin sizes used in Figure 8, substantially smaller than the 
typical DC signal. This quantity is independent of the primary charge Z. Since 
the DC signal increases as Z 2 , the night sky background contribution to the 
charge resolution AZ/Z decreases like 1/Z 2 . 
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Fig. 8. Direct and EAS Cerenkov light emitted from 100 TeV Z = 26 nucleus with 
a vertical trajectory. The Cerenkov light is observed at a mean radius of 80 m from 
the shower axis. Multiple scattering of the Cerenkov light due to Rayleigh and Mie 
scattering from a wind-driven aerosol model of the atmosphere is included. The 
vertical axis is the time delay of the arriving photons as discussed in the text. The 
horizontal axis is the arrival angle of the photons with respect to the vertical at the 
observing site. The scales below each panel give the photon intensities. 

4-1-2 Secondary Cerenkov Background 

Since the threshold energy for observing DC light increases linearly with in- 
creasing Z, the secondary Cerenkov light background, which is proportional to 
primary energy, also increases linearly. In principle, since the primary cosmic 
ray energy is measured from the secondary Cerenkov light, the amount of sec- 
ondary background light in the DC pixel bins could be estimated. This could 
be subtracted out to yield a pure DC light measurement. However, Poisson 
variations in the secondary Cerenkov background will generate fluctuations 
in the background-subtracted signal, thereby limiting the DC light measure- 
ment resolution. Consequently, the contribution of the secondary Cerenkov 
background to the charge resolution has the form 

AZ/Z oc ^ = Z^ 

In principle, optimal charge resolution is obtained by matching the detector 
pixel size and time resolution to the inherent width of the DC light emission. 
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In our simulation plots, we have used time bins of 500 ps and angular bins 
of 0.00375°, close to the optimal values. However, a conservative estimate of 
the magnitude of the secondary light yield has been determined from the 
simulations using a coarse 0.2° by 6 ns square integration window about the 
DC light yield position for a typical shower, in the absence of the DC light 
signal. This larger bin size is typical of existing state-of-the-art ground based 
Cerenkov observatories. 



4-1.3 Atmospheric Scattering 

The emitted DC light can possibly be scattered in the lower regions of the 
atmosphere. Mie and Rayleigh scattering are effects which scatter light from 
the EAS Cerenkov component into the DC light beam and also in to the 
gap in the time versus angle image (see Figure 3) which separates the EAS 
emission and the DC light. Rayleigh scattering by air molecules is essentially 
constant with respect to atmospheric conditions. The Rayleigh process scatters 
light to large angles (typically 90°) with respect to the original light beam 
direction, and therefore should provide very little effect on the separation 
between DC light and EAS Cerenkov light. Its main result is to increase the 
general background light level in every pixel rather than to favor a particular 
pixel. Less than 5% of the total Cerenkov light is scattered from the narrow 
image, essentially isotropically. The general increase in background light level 
due to Rayleigh scattering is negligible. 

Mie scattering is scattering from small suspended particles (aerosols) in the 
atmosphere, and is highly dependent on the aerosol composition, particle size 
distribution and height distribution. The Mie process acts to deflect the orig- 
inal photon beam out to angles of 5-10° off-axis, which is substantially more 
narrow than the characteristic Rayleigh angle. Because the Mie-scattered light 
is deflected to relatively narrow angles with respect to the original beam it 
represents a more serious potential to degrade the performance of the DC tech- 
nique. The amount of light scattered from the original beam depends strongly 
on the aerosol concentration. There are times when absolutely no light is Mie 
scattered ("molecular" atmosphere), and there are times when the atmosphere 
may be essentially opaque due to rain, dust storms, pollution, etc. If one is 
careful to select observation sites which have a large number of purely molec- 
ular atmospheric condition at night, one may be able to avoid this problem 
altogether. 

For example, the HiRes Fly's Eye observatory, located at the Dugway Proving 
Grounds, Utah, has a large fraction of its observing time (~ 90%) with viewing 
conditions similar to or better than the "standard desert atmosphere" [21]. In 
this wind-driven aerosol model, typically 5 — 20% of the total light beam is 
diverted uniformly out to angles of 5 — 10° off axis. Because the scattered 
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beam is wider than the original Cerenkov beam width, the general effect is 
to slightly elevate the background light levels in all pixels rather than in just 
the pixels relevant for the DC light measurement. Using small pixels, the 
effect of the Mie scattered light is substantially reduced. Figure 8 shows the 
level of background light, including atmospheric scattering, for a 100 TeV 
iron nucleus event. Clearly the DC emission remains prominent. At this level, 
the Mie scattering will likely affect the DC emission from the lighter nuclei 
(Z < 6). Heavier nuclei will be unaffected as their DC light is much stronger 
than the Mie scattered light from the EAS. 

For purposes of this paper, the resolution limitations due to atmospheric scat- 
tering of the secondary Cerenkov light is included in the secondary Cerenkov 
background calculation (Section 4.1.2). 

4-2 Detector Resolution Effects 

The detector-induced limitations to the charge resolution include the core 
position resolution, angular reconstruction resolution, and the error in the 
measured signal due to fluctuations in the DC signal photo-electron statistics 
(which depends upon the camera Quantum Efficiency (QE) and the Mirror 
area (A)). In the below calculations we assume A = 100 m 2 and QE = 25%. 

4-2.1 Angular Reconstruction Error 

For typical gamma-ray observatories like VERITAS [22] and HESS [23], a pri- 
mary trajectory angular reconstruction error of m 0.1° is expected. In these 
simulations, we assume the high-resolution cosmic ray detector has a similar 
angular reconstruction error. The charge resolution error is then computed by 
comparing average DC light yields at a mean radius of 80 m for an ensemble 
of vertical simulated showers with another ensemble of showers with identical 
charge and energy, but with a primary zenith angle trajectory of 0.1°. The 
resulting error is independent of Z, and very small (AZ/Z xs 0.34%). 

4-2.2 Gore Position Error 

Given a gamma-ray observatory like VERITAS [22] with an array of 10 m diam- 
eter primary mirrors, the core position can only be localized to approximately 
5 meters. Assuming a 5 meter core resolution error, we examine the change 
in the DC light yield as a function of distance from the shower core using the 
Monte Carlo simulation. Figure 7 shows the results of simulations of several 
thousand vertically incident 100 TeV iron nuclei viewed in a 0.2° angular by 6 
ns delay-time bin centered on the DC emission region. The dashed histogram 
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shows the DC light observed in an annulus of 67 — 94 m (mean radius 80 m) of 
the particle path, the solid histogram shows the DC light observed within an 
annulus of 47 — 67 m (mean radius 57 m). Clearly both of these observations 
have well defined peaks. Simulations of different Z nuclei show virtually the 
same ratio of the peak DC light yield at mean radii of 57 m to 80 m. Using 
a linear interpolation of the light yield variation with distance to the shower 
core, one derives a charge resolution error of 3.1% due to core position un- 
certainty. This charge resolution error is independent of the primary charge 
Z. 



4-2.3 Photon Statistics 

Fluctuations in the photon statistics play an important role for low DC light 
emission levels (at small Z). Since the signal is proportional to Z 2 , the charge 
resolution scales like 

AZ/Z oc ^ = 1/Z. 

The magnitude of signal is determined from the Monte Carlo Simulation, 
assuming the above specified values for QE and A. 



4-3 Hadronic Interaction Fluctuations 



Because the primary particle will eventually suffer hadronic interaction, the 
length of the DC light emitting region changes from event to event with the 
fluctuation in the depth of the first interaction. This leads to a subsequent 
fluctuation in the DC light yield. If the interaction point can be identified in 
the track image (appearing as a sudden drop in emission rate in a sufficiently 
finely pixelation detector), the integrated yield could be normalized to path 
length. However, such an analysis in reality is complicated by a finite detector 
resolution, photon statistics, and the fact that the nucleus does not usually 
suffer a catastrophic collision in the first interaction. Typically, a nucleus of 
charge Z will break up into a charge Z — 2 nucleus and a Helium (Z = 2) 
nucleus, and therefore continue to emit DC light at only a modestly reduced 
rate until it suffers further collisions. Fluctuations in this fragmentation pro- 
cess therefore affect the integrated DC light yield in a fashion for which there 
may be no simple correction available. 

The fluctuation effects are best illustrated by examining characteristics of a 
large ensemble of simulated showers. A scatter plot (Figure 9 ) illustrates the 
distribution of integrated DC light measured between 67 — 94 m from the 
shower core as a function of the depth of the first nuclear interaction of the 
primary. The scatter plot has essentially three regions of interest, depending 
upon the depth of the first interaction. 
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For cosmic rays which interact less than 7 g/cm 2 deep (i.e., at altitudes above 
~ 40 km), there are large fluctuations in the DC light yield. Even nuclei 
with the same first interaction depth suffer from large fluctuations. The DC 
light yield fluctuation is due mainly to the details on how the first interaction 
proceeds. If the initial interaction yields a large number of low Z fragments, 
the light yield is low. If the interaction yields only 1 or 2 low Z fragments and 
the original primary only loses a small amount of charge, then the DC light 
yield is high. This region (interaction depth < 7 g/cm 2 ) may be referred to as 
the 'fragmentation region' as the magnitude of the DC light yield in this region 
is dominated by the details of the fragmentation process. The fragmentation 
region is illustrated in lower region of Figure 9. 

For cosmic rays which interact between 7—18 g/cm 2 in the atmosphere, there 
is a tight linear relationship between depth of first interaction and the DC light 
yield. In essence, as the particle penetrates deeper into the atmosphere, the 
total amount of DC light emitted by the full nucleus begins to dominate any 
DC light fluctuations generated by fluctuations in the fragmentation process. 
The DC light yield can be though of purely geometrically, as the length of 
the path traveled between 7—18 g/cm 2 before the first interaction in this 
region. This region may be called the 'cross-section' region as the DC light 
yield in this region is dominated by the value of the inelastic cross-section. 
The cross-section region is illustrated in middle region of Figure 9. 

For first interactions occurring deeper than 18 g/cm 2 (i.e., altitudes below 
~ 30 km), the light yield is essentially independent of the interaction depth. 
This phenomena arises because the DC light is emitted at a specific emission 
angle. The emission angle, when coupled with the height of emission, defines 
the radius at which the DC light will be observed. When the particle passes 
beyond 18 g/cm 2 without interacting, all the DC light emitted beyond this 
depth falls outside the 67—94 m radial bin used in this plot. Consequently, once 
the nucleus has passed through 18 g/cm 2 without interacting, all light that 
can be emitted into the 67 — 94 m annulus has been emitted, independent of 
where the subsequent first interaction takes place. This region may be referred 
to as the 'saturation region', as the DC light yield is essentially saturated to 
its maximum value. The saturation region is illustrated in upper region of 
Figure 9. 

If a detection annulus closer to the particle path is used, the DC collected will 
be emitted higher in the atmosphere, as shown in Figure 2. The altitude or 
grammage boundaries between the fragmentation, cross-section, and satura- 
tion regions are therefore functions of the radial distance to the shower core. 
In this case a particle can interact higher in the atmosphere (still below 40 
km), existing in the cross-section region for the outer radius (80 m) but still 
be in the saturation region for this inner radius (57 m). In addition, the overall 
light yield in these regions decrease as one decreases the observation radius. 
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Fig. 9. Correlation of DC emission with pathlength in atmosphere before first in- 
teraction for 100 TeV Iron primary [Z = 26). Horizontal Axis: DC light intensity 
at 67 — 94 m (mean radius 80 m) from shower core (photons/m 2 ). Vertical Axis: 
Pathlength in atmosphere before first interaction (g/cm 2 ). 



Since particles can exist in the saturation region for one radius and be in the 
cross-section region for a larger radius, it is possible to extract additional in- 
formation by examining the ratio of DC light emitted in different radial bins. 
For instance, for a particle in the saturation region at 57 m but in the cross- 
section region at 80 m, the ratio of DC light measured at the 80 m distance 
to the DC light measured at the 57 m distance is strongly correlated with the 
depth of the first interaction. Consequently, a model-independent measure- 
ment of the inelastic nuclear cross-section is possible, even with moderately 
coarse camera pixels. Determination of size of the fluctuations in the frag- 
mentation process may also be possible using three radial measurements: an 
outer radius measurement in the fragmentation region, a middle radius mea- 
surement in the cross-section region, and an inner radius measurement in the 
saturation region. 

An important concept is that all possible interaction fluctuations can only 
result in a decrease in the observed DC light level. This makes DC light dis- 
tribution in Figure 10 asymmetric at the lower DC light intensity level. For 
example, although Z = 64 nucleus might occasionally yield a very low light 
DC signal comparable to a Z = 26 or Z = 40 nucleus, a smaller charge (e.g., 
Z = 40) nucleus cannot emit the DC light intensity of a Z = 64 nucleus. 
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Fig. 10. Distribution of DC light emitted at mean radius 80 m from shower core for 
several thousand simulated showers at each primary Z. Horizontal Axis: DC light 
intensity (photons/m 2 ). Vertical Axis: Probability (1.0 =100%) 



Consequently interaction fluctuations cannot result in an overestimation of 
the particle charge, it can only underestimate it. 

The interaction fluctuation contribution to the overall resolution is determined 
from the width of each narrow peak distribution in Figure 10. The resulting Z 
dependence of the charge resolution is found to be well described by a power 
law 

AZ/Z oc Z-°- 73 . 



4-4 Overall Resolution 



The charge resolution expected from the DC technique is limited by the com- 
bination of the various effects described above. In order to be conservative, we 
have assumed a detection scheme with an effective collecting area of 100 m 2 , 
a core location capability of 5 m, a time resolution of 6 ns, and an angular 
pixel size of 0.2°. Figure 11 shows the charge resolution expected resulting 
from these effects as a function of charge Z. For low charges the resolution is 
dominated by secondary Cerenkov light from the EAS. At higher Z the core 
resolution provides the charge resolution limitation. The overall resolution is 
calculated to be AZ/Z ~ 5% for Z > 10, essentially independent of charge. 
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Fig. 11. The expected charge resolution AZ/Z for a detector of effective area 100 
m 2 and core position resolution 5 m. Horizontal Axis: Primary Charge Z. Vertical 
Axis: Charge Resolution AZ/Z{%) 

In principle, improved charge resolution could be obtained by matching the 
detector pixel size and time resolution to the inherent width of the DC light 
emission (about 500 ps timing and 0.00375° wide pixels). Since the resolution 
function for high Z is dominated by core position error contribution, this would 
not improve the resolution at large Z. However, it could substantially improve 
the resolution for Z < 10 as it would reduce the effects of the secondary 
Cerenkov light in the charge resolution. 



5 Detector Design and Applications 

A possible design for a DC imaging observatory can be constructed from the 
simulation work results. Such an observatory comprises a number of large 
(>10 m diameter) front surface reflecting dishes in fixed mounts. The mounts 
point vertically for low energy composition measurements, but are adjustable 
to view at large zenith angles to increase detection area for measuring the 
higher energy cosmic ray heavy nuclei. The dishes are separated by a distance 
of ~ 80 m for sea level operation, placed on a hexagonal grid pattern. 

Each dish has an isochronous surface in order to preserve the incoming wave- 
front timing to approximately 2 ns or better. This timing requirement provides 
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additional separation between the DC light and the secondary Cerenkov light 
as shown in Figure 3. A highly pixelated optical camera built from small, fast 
photomultiplier tubes is used at the focal plane to image the Cerenkov light 
from the cosmic ray shower. In order to have the best resolution and signal- 
to-noise ratio for the DC light, small pixels (typically 0.01°) are preferred. 
However, the DC light measurements can be made with large pixels, up to 
0.1 — 0.15°. A field of view of 5 — 10 degrees is needed to provide a sufficient 
collection aperture. For a 10° x 10° field of view, and a fiducial area defined 
between 50 — 115 m of a detector, each dish has an effective collection aperture 
of ~ 1000 m 2 sr. 

Light striking the individual pixels is read out by a fast FADC system, with 
500 MHz sampling rate or faster, enabling an image to be constructed in time 
versus/angle like those shown in Figure 3. A key advantage of the DC method 
over the observations of Cerenkov light high in the atmosphere is that the 
triggering of the system is relatively simple. Since a significant part of the 
EAS will develop in the field of view of the telescopes, this large light sig- 
nal can be used to provide a simple trigger scheme. A Cerenkov light signal 
produced by an EAS with energies > 10 TeV can be reliably discriminated 
from the night sky by simple logic on the camera pixels for a 10 m mirror 
size. Using individual-tube constant fraction discriminators, hit patterns from 
multiple pixels are combined using fast logic tables to look for sequences of 
pixels with characteristics similar to a Cerenkov image. Once a mirror is trig- 
gered, it retrieves a history of photon time slices from the FADC system for 
about 20 — 50 nanoseconds around the trigger time to look for the delayed 
DC light signal. Neighboring telescopes in the array would also be read out 
to look for coincidence measurement of the same Cerenkov image from dif- 
ferent observation angles. The data from multiple telescopes is combined to 
determine the shower geometry, energy, and DC light content. Importantly, 
events can be required to contain consistent amounts of DC light intensity and 
location in multiple telescopes. This procedure can reject essentially all light 
signal contamination from local sources, for example local muons, and verify 
the level of DC emission from independent measurements. 

The total detection aperture is not only affected by the physical size of the 
obervatory and the field of view of each camera; it is also affected by addi- 
tional constraints needed to observd the DC, light such as requiring the shower 
core to strike in an annulus between 50-115 m from the individual detectors. 
Because we are using an array of multiple telescopes, in general the shower 
core will fall within an acceptable distance from more than one telescope over 
most of the area surrounding the array. Additional considerations, such as field 
of view, also restricts the detection aperture. These considerations are simi- 
lar to those imposed on typical ground-based gamma ray observatories (e.g. 
VERITAS), and consequently we expect a similar detection aperture to these 
devices (i.e. several thousand m 2 sr). However, an accurate determination of 
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the detector acceptance, including energy and charge dependence of the aper- 
ture, can only be accomplished through detailed Monte Carlo simulation. This 
simulation must include a detailed model of an assumed detector configura- 
tion (array spacing, mirror area, pixel size, electronic triggering charateristics, 
etc). Monte Carlo calculations for a typical strawman observatory configu- 
ration, modeled after the VERITAS observatory, are underway and will be 
presented in a subsequent paper. 

We note that the next generation of imaging Cerenkov telescope arrays, in- 
cluding VERITAS [22], HESS [23], and CANGAROO[24] all begin to approach 
the strawman design, and will likely be able to provide a first measurement 
of DC light. An intriguing possibility is to use the DC light to reject back- 
ground in these experiments. The ability to distinguish between the Cerenkov 
signals produced by electrons and gamma rays in the 10 — 500 GeV energy 
range is important for substantial improvement in the sensitivity of ground- 
based gamma-ray astronomy. In this energy range, protons and heavy nuclei 
do not produce secondary particles with sufficient energies to generate sub- 
stantial Cerenkov light. The majority of the background events for gamma ray 
instruments in this region are produced by cosmic ray electrons, which also 
generate pure electromagnetic air showers, with EAS Cerenkov image charac- 
teristics identical to a gamma ray primaries. The capability to identify even 
a small fraction of the electron events could make significant improvements 
in the sensitivity of gamma-ray telescopes in this energy range. The upper 
left hand panel of Figure 6 shows an angle/time image of a simulated vertical 
gamma-ray shower at 10 TeV. The upper left panel shows a vertical electron 
at 100 GeV. This second event could be identified as an electron if the faint 
DC image arc at the left of the panel is detected. The intensity in this arc is 
~ 5 times the expected night sky background. 

A similar background occurs in the 500 GeV —10 TeV range for protons which 
transfer most of their primary energy to one or two ir° secondaries in their 
first interaction, leaving little energy for further hadronic interaction (i.e., 
generation of which decay to muons). These types of events represents 
a substantial amount of the residual cosmic ray background after applying 
standard image analysis cuts to select gamma rays. Since these proton-induced 
events are pre-selected to have image characteristics similar to gamma rays, 
additional parameters must be used to reject these events. The lower left 
hand panel in Figure 6 shows an event in which the DC signal from a proton 
has transferred most of its energy to an initial n° during the interaction. In 
comparison with the gamma-ray events these type of events might possibly be 
rejected using finely-grained timing and angular imaging. 
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6 Conclusions 



We have discussed a new experimental technique which can potentially yield 
excellent charge resolution measurements (AZ/Z < 5% for Z = 26) for 
ground-based observations of high energy cosmic rays. The technique relies 
upon the observation of the direct Cerenkov light emitted by the primary nu- 
cleus before the first nuclear interaction with the Earth's atmosphere. The 
experimental technique works over an energy range in the TeV-PeV energy 
range, with the width of the energy window growing like Z for heavy nuclei. 
A dedicated observatory could be built to observe cosmic rays using his tech- 
nique, while the next-generation ground-based imaging Cerenkov telescopes 
such as VERITAS, HESS, and CANGAROO will begin to approach the sen- 
sitivity required to begin initial observations. 

The average intensity and fluctuations in the DC light yield have been ex- 
amined using Monte Carlo simulations. The yield contains the expected Z 2 
dependence with fluctuations that are easily understandable in terms of the 
radial distance of observation, variations in the depth of first hadronic inter- 
action, and the details of the subsequent fragmentation process. Indeed, DC 
light measurement appears to provide a method to measure the primary charge 
independent of any hadronic interaction or fragmentation model. This could 
potentially solve the most difficult systematic problem in the ground-based 
studies of TeV/PeV cosmic rays: the difficulty in distinguishing between ac- 
tual changes in composition and systematic shifts of nuclear interaction models 
with energy. The high resolution nature of this measurement might allow the 
detection of changes in the propagation pathlength distribution at high energy. 
The pathlength distribution at PeV energies is unknown, but has a extremely 
strong effect on the predicted composition at the knee of the all-particle spec- 
trum [25]. 

The DC light measurement technique could provide researchers with a "tagged 
particle beam" of known energy and composition. By examining DC light yield 
at various radial bins, one may be able to extract nucleon-air cross-section in- 
formation as a function of primary mass and energy, independent of an inter- 
action model. Fluctuations in the DC light yield for small interaction depths 
could provide information concerning the fluctuations in the fragmentation 
processes of nucleus-air interactions. 

The DC light measurement technique may provide additional sensitivity to 
ground-based imaging air Cerenkov gamma-ray detector for the rejection of the 
dominant electron background in the 10 — 500 GeV energy range, and residual 
proton background in the 1 — 10 TeV energy range. Such an improvement 
could potentially increase the point source sensitivity of these telescopes by 
an order of magnitude or more. 
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The unique properties of DC light (versus EAS light) may allow the observa- 
tion of exotic cosmic ray particles, such as "strange quark matter" (see, e.g., 
[26,27]) or magnetic monopoles [28]. These particles, which are thought to 
have large effective charges, would provide an extremely strong DC signature 
due to the Z 2 dependence of the light yield. 

A technique which combines the collection power and logistical advantages of 
ground-based cosmic ray detectors with the high-precision charge resolution 
of balloon or satellite-borne experiments would be a valuable asset in the field 
of high-energy astrophysics. Such a technique, if scientifically and financially 
viable, has the potential to make great strides in the determination of cosmic 
ray composition across the "knee", and subsequently, foster advances in the 
theory of cosmic ray origins. We believe that the development of imaging at- 
mospheric Cerenkov telescopes has provided the means to target the Cerenkov 
light emitted directly from primary nuclei prior to their hadronic interaction 
in the atmosphere. A direct Cerenkov experiment, implemented as an array 
of 10 m reflectors imaged by high-resolution cameras with fast ADC systems 
may prove to be the next step towards understanding the nearly century-old 
mystery of cosmic rays. 
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